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Abstract 

We study capability of f(R) gravity models to allow crossing the phantom boundary 
in both Jordan and Einstein conformal frames. In Einstein frame, these models are 
equivalent to Einstein gravity together with a scalar field minimally coupled to gravity. 
This scalar degree of freedom appears as a quintessence field with a coupling with the 
matter sector. We investigate evolution of the equation of sate parameter for some 
cosmologically viable f(R) gravity models in both conformal frames. This investigation 
(beyond mere theoretical arguments) acts as an operational tool to distinguish physical 
status of the two conformal frames. It shows that the two conformal frames have not the 
same physical status. 

PACS Numbers: 98.80.-k 



1 Introduction 

For explaining accelerated expansion of the universe, there is a class of models in which one 
modifies the laws of gravity whereby a late-time acceleration is produced. A family of these 
modified gravity models is obtained by replacing the Ricci scalar R in the usual Einstein- 
Hilbert Lagrangian density for some functions f(R) [1] [2]. It is well-known that [3] [4] the 
gravitational field equations derived from such fourth order gravity models can be conformally 
transformed to an Einstein frame representation with an extra scalar degree of freedom. This 
scalar degree of freedom may be regarded as a manifestation of the additional degree of freedom 
due to the higher order of the field equations in the Jordan frame. This feature opens some 
questions about a possible correspondence between this geometric scalar field and the usual 
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scalar fields such as quintessence and phantom fields. Along this line, we intend in the present 
contribution to explore phantom behavior of some f(R) gravity models in both Jordan and 
Einstein conformal frames. In particular, we shall show that behaviors of the scalar field 
attributed to f(R) models are similar to those of a coupled quintessence, a quintessence which 
interacts with the matter sector. 

The plan of our paper is as follows : In section 2, we consider the cosmology of a scalar 
field with minimal coupling to gravity. We review the observation that a quintessence field 
with a minimal coupling can not lead to crossing the phantom boundary. In section 3, we 
shall consider f(R) gravity in both Jordan and Einstein conformal frames. We show that the 
generic feature of these models is that the scalar partner of the metric tensor in the Einstein 
frame is similar to a quintessence rather than a phantom field. The interaction of this scalar 
degree of freedom with the matter sector plays a key role for crossing the phantom barrier [5] . 
We then compare phantom behavior of some viable f(R) models in the two conformal frames. 
Our results are summarized in section 4. 



2 Minimally coupled Scalar Field 

The simplest class of models that provides a redshift dependent equation of state parameter 
is a scalar field ip minimally coupled to gravity whose dynamics is determined by a properly 
chosen potential function V(ip). Such models are described by the action t 

S v = I J d 4 xV=9 {\R - a g^d^ip - 2V(<p)) + S m (g^, (1) 

where k = 8nG with G being the gravitational constant, g is the determinant of g^ u and R 
is the curvature scalar. Here S rn is the action of dark matter which depends on the metric 
g^ v and some dark matter fields ip. The constant a can take a = +1, —1 which correspond to 
quintessence and phantom fields, respectively. The distinguished feature of the phantom field 
is that its kinetic term enters (1) with opposite sign in contrast to the quintessence or ordinary 
matter. The action (1) gives the Einstein field equations 

GV = k(T£, + Tp (2) 

with | 

T ^v = « V M y?V^ - -a ^V 7 ^V°V - 9^V((p) (3) 

Here T™ is the stress-tensor of the matter system defined by 



^ ,/=g Sg^ 



The two stress-tensors and T™ are separately conserved 

V"7£ = V%™ = (5) 



1"We work in units in which % = c = 1 and the signature is (—,+,+,+). 
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We assume a spatially flat homogeneous and isotropic cosmology described by Friedmann- 
Robertson- Walker (FRW) spacetime 

ds 2 = -dt 2 + a 2 (t)(dx 2 + dy 2 + dz 2 ) (6) 

where a(t) is the scale factor. In this cosmology, the gravitational field equations (2) become 

3H 2 = k(p^ + p m ) (7) 

2H = -k[(u v + l)p v + p m ] (8) 
where H = ^ is the Hubble parameter and 

Pp = + V (v) ' Pf = ~ V ^ ^ 



\aup 2 - V(tp) 
\aip 2 + V(ip) 



= C 2 , ; : do) 



The conservation equations (5) take the form 

p m + 3Hp m = (11) 

p\ + 3H(w v + l)p v = (12) 

In the case of a quintessence field (a = +1) with V((p) > the equation of state parameter 
remains in the range — 1 < u v < 1. This is also true for a phantom field (a = —1) with a 
negative potential V(<p) < 0. In the limit of small kinetic term (slow-roll potentials [6]), it 
approaches = — 1 but does not cross this line. For a — +1, we have 

+ P v = (u v + l)Pp = cup 2 > (13) 

The phantom barrier can be crossed by a phantom field (a < 0) with V(ip) > when we have 
2|V(<£>)| > (p 2 . This situation corresponds to 

Pp > , Pv < , V(<p) > (14) 

In this case, (ou^ + l)p v < and for a sufficiently negative pressure p v the equation (8) gives 
H > 0. Let us look at this situation in a different way. We combine (8) with (11) and (12) 
which leads to 

2H = ^ ! {Pp + p m ) (15) 

In an expanding universe p m < 0. In the case of a quintessence (uu^ + > and then the 
equation (12) implies that p\ < which leads to H < 0. On the other hand, for a phantom 
field we have a < and (u;^ + l)p v < 0. In this case, p^ > and for an appropriate potential 
function V(tp) the first term on the right hand side of the equation (15) dominates and the 
equation then results in H > 0, crossing the phantom barrier. The merit of (15) is that it 
implicitly shows that there is a possibility for an interacting quintessence field to cross the 
phantom boundary. We will return to this issue in the next section. 

Here it is assumed that the scalar field has a canonical kinetic term ±|<p 2 - It is shown [7] that 
any minimally coupled scalar field with a generalized kinetic term (k-essence Lagrangian [8]) 
can not lead to crossing the PDL through a stable trajectory. However, there are models that 
employ Lagrangians containing multiple fields [9] or scalar fields with non-minimal coupling [10] 
which in principle can achieve crossing the barrier. 
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3 f(R) Gravity 



The f(R) gravity models are based on a Lagrangian density which depends in a nonlinear way 
on the curvature scalar. In these models the dynamical variable of the vacuum sector is the 
metric tensor and the corresponding field equations are fourth order. This dynamical variable 
can be replaced by a new pair which consists of a conformally rescaled metric and a scalar 
partner. Moreover, in terms of the new set of variables the field equations are those of General 
Relativity. The original set of variables is commonly called Jordan conformal frame and the 
transformed set whose dynamics is described by Einstein field equations is called Einstein 
conformal frame. 

In general, the mathematical equivalence of these two conformal frames does not imply a 
physical equivalence. The physical status of these conformal frames is an open question which 
has not been completely solved yet. In fact, there are different categories according to various 
authors' attitude towards the issue of which frame is the physical one [4] [11]. Based on the 
interaction of the scalar degree of freedom with the matter sector, some authors argue that the 
two conformal frames are physically equivalent provided that one accepts the idea that the units 
of mass, length and time are varying in the Einstein frame [12]. According to this idea, physics 
must be conformally invariant and symmetry group of gravitational theories should be enlarged 
to include not only the group of diffeomorphisms but also conformal transformations. There 
are also arguments against this idea and in favor of one of the conformal frames. In particular, 
some authors go beyond mere theoretical arguments and take an appropriate phenomenon to 
determine if the two frames are physically equivalent [13]. 

One important issue in the context of cosmological viability of an f(R) gravity model is 
capability of the model to exhibit phantom behavior due to curvature corrections. We shall 
explore this issue in both Jordan and Einstein conformal frames and take it as a relevant issue 
to distinguish the physical status of the two conformal frames. 

3.1 Jordan Frame Representation 

The action for an f(R) gravity theory in the Jordan frame is given by 



Stability issues should be considered to make sure that an f(R) model is viable [14]. In partic- 
ular, stability in matter sector (the Dolgov-Kawasaki instability [15]) imposes some conditions 
on the functional form of f(R) models. These conditions require that the first and the second 
derivatives of f(R) functions with respect to the Ricci scalar R should be positive definite. 
The positivity of the first derivative ensures that the scalar degree of freedom is not tachyonic 
and positivity of the second derivative tells us that graviton is not a ghost. 
The field equations can be derived by varying the action with respect to the metric tensor 




(16) 



f{R)R, v - \g^f{R) - V,V u f(R) + g^nf(R) = kT; 



(17) 



where prime denotes differentiation with respect to R. The trace of (17) is 

f'(R)R - 2f(R) + 3af(R) = kT m 



(18) 
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where T m = g^T™. It should be noted that when f(R) = R, the equation (18) reduces to 
R = —kT m which is the corresponding trace equation in General Relativity. In f(R) gravity, 
□/'(i?) does not vanish and contrary to General Relativity the Ricci scalar relates differentially 
to the trace of the matter system. This is an indication of the fact that one finds a larger 
variety of solutions in f(R) modified gravity models. 
The field equations can also be written in the following form 



= K ff (T™ + T™) (19) 

where fc e // = k/f'(R) and 



l r l 



tst = i^(f(R) - Rf(R))9r + v M v„/'(i2) - g^rm m 

We apply the field equations to a spatially flat FRW spacetime described by (6). We assume 
that the matter system is described by a pressureless perfect fluid with energy density p m . The 
field equations become 

3H 2 = k eff ( Pm + p DE ) (21) 
2H + 3H 2 = -k eff p DE (22) 

where 

Pde = ~ f(R)) ~ 3HRf"(R)] (23) 

pde = \\\U\R) - Rf'(R)) + Rf"(R) + R 2 f"'(R) + 2HRf"(R)] (24) 

Then the effective equation of state parameter is 

^// = - P Z = ^A7^ [kf(R) - Rf'(R)) + Rf(R) + R 2 f"(R) + 2HRf{R)\ (25) 
Pm + Pde oH z / (K) I 

where we have used (21). This expression allows us to find those f(R) functions that fulfill 
uj e ff < —1. In general, to find such f(R) gravity models one may start with a particular f(R) 
function in the action (16) and solve the corresponding field equations for finding the form of 
H(z). One can then use this function in (25) to obtain uu e ff(z). However, this approach is 
not efficient in view of complexity of the field equations. An alternative approach is to start 
from the best fit parametrization H(z) obtained directly from data and use this H(z) for a 
particular f(R) function in (25) to find u> e ff(z). In this approach, one needs a consistency 
check to ensure that the given parametrization is compatible with the f(R) function. We will 
follow the latter approach to find f(R) models that provide crossing the phantom barrier. We 
begin with the Hubble parameter H which its derivative with respect to cosmic time t is 

H=-~ & (26) 
a a 

Combining this with the definition of the deceleration parameter 

<m = ~ (27) 
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gives 

H = -{q + l)H 2 (28) 

One may use z = — 1 with z being the redshift, and the relation (27) to write (28) in its 
integration form 

H(z) = H exp[ f (1 + q(u))d\n(l + u)} (29) 
Jo 

where the subscript "0" indicates the present value of a quantity. Now if a function q(z) is 
given, then we can find evolution of the Hubble parameter. Here we use a two-parametric 
reconstruction function characterizing q(z) [16] [18], 

1 q 1 z + q 2 ( . 

q{z) = 2 + JTTW m 

One of the advantages of this parametrization is that q(z) — > \ when z ^> 1 which is consistent 
with observations. Moreover, The behavior of q(z) in this parametrization is quite general. If 
qi > and q 2 > 0, then there is no acceleration at all. It is also possible that the Universe 
has been accelerating since some time in the past. If q± < and q 2 > —1/2, then it is possible 
that the Universe is decelerating and has past acceleration and deceleration. In other words, 
it has the same behavior as the simple three-epoch model used in [16] [19]. The values of qi 
and q 2 and the behavior of q(z) can be obtained by fitting the model to the observational data. 
Fitting this model to the Gold data set gives q x = lA7t\H an d <?2 = -1.46 ±0.43* [18]. Using 
this in (29) yields 

H(z) = H (l + zf> exp[| + |^|] (31) 

For the metric (6), we have R = 6(H + 2H 2 ) and therefore R = 6(H + AHH). In terms of the 
deceleration parameter, one can write 

R = 6(1 - q)H 2 (32) 

and 

R = 6H 3 {2(q 2 - 1) - |} (33) 

R = QH 3 {Q(q 2 - 1)| + 4qq - 2^ - |} (34) 
The latter two equations are equivalent to 

R = 6H 3 {2(q 2 - 1) + (1 + zA (35) 

az 

R = -6(z + l)H 3 {3^[2(q 2 - 1) + (z + l)g] + H[(4q + l)g + (z + 1)0]} (36) 
From the equations (30) and (31) we can also write 

dq (<?i - 2q 2 ) - q x z 



dz {z+lf 



(37) 



■l-Some other recent samples may be found in [17]. 
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fq _ 2 qi (z-2) + 6q 2 

= Fo(1 + z) [ 2 + (^Ti7 ] exp[ 2 + 2(7TTF ] (39) 

It is now possible to find R, R and i? in terms of the redshift, and then for a given f(R) 
function, the relation (25) determines the evolution of the equation of state parameter u e ff(z). 
As an illustration we apply this procedure to some f(R) functions. Let us first consider the 
model [20] [21] 

f(R) = R + XRo& n (40) 

Here Rq is taken to be of the order of Hq and A, n are constant parameters. In terms of the 
values attributed to these parameters, the model (40) is divided in three cases [21]. Firstly, 
when n > 1 there is a stable matter-dominated era which does not follow by an asymptoti- 
cally accelerated regime. In this case, n = 2 corresponds to Starobinsky's inflation and the 
accelerated phase exists in the asymptotic past rather than in the future. Secondly, when 
< n < 1 there is a stable matter-dominated era followed by an accelerated phase only for 
A < 0. Finally, in the case that n < there is no accelerated and matter-dominated phases 
for A > and A < 0, respectively. Thus the model (40) is cosmologically viable in the regions 
of the parameters space which is given by A < and < n < 1. 
For the model (40), one can write 

f\R) = l + nA(^ r i 

f»(R) = n(n - l)\H Q \^Y-* (41) 

f"\R)=n{n-l){n-2)\H^\^r^ 

Putting these results together with (31), (32), (35) and (36) into (25) gives u e ff(z). Due to 
complexity of the resulting u> e ff(z) function, we do not explicitly write it here and only plot it 
in fig. la for some parameters. The figure indicates that u e ff(z) crosses the phantom boundary 
for some values of the parameters. 

Now we consider the model presented by Starobinsky [22] 

f(R) = R-lR {l-[l + Ar m } (42) 

where 7, m are positive constants and Rq is again of the order of the presently observed 
effective cosmological constant. For this model, one obtains 



R R 
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f(R) = l-2m 1 ^-Jl + ^ I r m - 1 



H 2K R 4 - 



f"(R) = 2m 7 iJ - 2 (l + ^)-™- 1 [2(m + 1)^(1 + j*)" 1 - 1] (43) 
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f"{R) = 4m(m + 1) 7 ^(1 + 3 [3 - (2m + 1)^] 

The corresponding uj e ff(z) function is plotted in fig. 2a. The figure indicates that the Starobin- 
sky's model also realizes crossing the phantom boundary. 

3.2 Einstein Frame Representation 

Any f(R) gravity model described by the action (16) can be recast by a new set of variables 

9iiv = ft 9txv (44) 

= ^^lnft (45) 

where Q = f'(R) and /3 = ^/|. This is indeed a conformal transformation which transforms 
the action (16) to the Einstein frame [3] [4] 

Sef = \I { l R - dA - + S ™(9^ 2 ^, ■*!>) (46) 

In the Einstein frame, is a minimally coupled scalar field with a self-interacting potential 
which is given by 

V(m) = « (47) 

Note that the conformal transformation induces the coupling of the scalar field <fi with the 
matter sector. The strength of this coupling (3, is fixed to be and is the same for all types 
of matter fields. 

Variation of the action (46) with respect to g^ u and gives, 



GV = KT*, + fp (48) 



U<f>-^- = -/3Vk f m (49) 



where 

Ttu = V M V v <t> - i^,V 7 - V{4>)g^ (51) 

Here T™ and T^ u are stress-tensors of the matter system and the minimally coupled scalar 
field 0, respectively. The trace of (48) is 

V 7 0V 7 + 4V((f>) - R/k = f m (52) 

which differentially relates the trace of the matter stress-tensor T m = g^f™ to R. It is 
important to note that the two stress-tensors T™ and T$ v are not separately conserved. Instead 
they satisfy the following equations 

V"7jJ = -V%t = f3 Vk V„0 f m (53) 
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We apply the field equations to a spatially flat FRW metric which in the Einstein frame is 
given by 

ds 2 = -dP + a 2 (t)(dx 2 + dy 2 + dz 2 ) (54) 

where a = Vfia and dt = Vt^dt. To do this, we take T™ and T^ u as the stress-tensors of 
a pressureless perfect fluid with energy density p m , and a perfect fluid with energy density 
P<j> — \ ( t )2 + y(0) an d pressure = |0 2 — V((f>), respectively. In this case, for the metric (54) 
the equations (48) and (49) take the form § 

3H 2 = k(p^ + Pm ) (55) 

2H = -fc[(a^ + l) P4> + p m ] (56) 

+ 3tf + = -pVk Pm (57) 

d(p 

The trace equation (52) and the conservation equations (53) give, respectively, 

02 + R/k _ W ^ = pm (58) 

p m + 3H Pm = Q (59) 

p + 3#(w + 1)P0 = ~Q (60) 

where 

Q = {3Vk<P Pm (61) 

is the interaction term. This term vanishes only for = const., which due to (45) corresponds 
to the case that f(R) linearly depends on R. The direction of energy transfer depends on the 
sign of Q or 0. For > 0, the energy transfer is from dark energy to dark matter and for 
< the reverse is true. 

The formulation of an f(R) theory in the Einstein frame indicates that such a modified gravity 
model has one extra scalar degree of freedom compared with General Relativity. This suggests 
that it is this scalar degree of freedom which drives late-time acceleration or crossing the 
PDL in a cosmologically viable f(R) gravity. This opens questions about the role of this 
scalar degree of freedom as quintessence or phantom. Comparing the sign of the kinetic terms 
of the scalar field in the actions (46) and (1) immediately reveals that this scalar degree of 
freedom actually appears as a quintessence rather than a phantom field. It is well-known that 
a minimally coupled quintessence field may lead to an accelerated expansion but it can not 
lead to crossing the PDL. This may lead one to conclude that crossing the PDL can not take 
place in models such as (46). However, there is a basic difference between the two actions 
(1) and (46). In the former there is no interaction between the two fluids while in the latter 
the scalar field interacts with the matter sector. In the following we shall show that this 
interaction is actually responsible for a possible crossing of the PDL in a viable f(R) model. 
To do this, we first combine (59) and (60) with (56) to obtain 

2H = ^(p<p + Pm) (62) 



3 Hereafter we will use unbarred characters in the Einstein frame. 
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Apart from some similarities between the latter and the equation (15), there is an important 
difference between the two equations concerning the interaction between the scalar field <fi and 
the matter part. In general, we can consider two different cases ^ : 

a) Firstly, Q > which corresponds to energy transfer from dark energy (or the scalar degree 
of freedom </>) to (dark) matter. Due to the fact that appears in our case as a quintessence, 
we have always p^oo^ + 1) > for both signs of Q. Thus in the case that Q > 0, p^ is always 
a decreasing function in an expanding universe (H > 0). From (59) one infers that p m can 
take both positive and negative signs depending on the relative magnitudes of Q and 3Hp m . 
The equation (62) implies that when p m > one can consider the possibility that H > and 
crossing the phantom barrier. 

b) Secondly, Q < corresponds to a reversed sign of energy transfer. In this case p m is definitely 
negative in an expanding universe. Instead p^ can be negative and positive depending on the 
relative magnitudes of p^ioj^ + 1) and Q. The case that p^ > may lead the accelerating 
expansion to cross the PDL. 

This argument emphasizes the role of interaction in the sign of H even though it does not 
appear explicitly in the equation (62). The effective equation of state parameter is defined by 

u eff = = -^+- (63) 

Peff P<f> + Pm 

We may write it in a more appropriate form. To do this, we can use (55), (59) and (60) that 
leads to 

k k 3 

^e// = -1 - ^jp(f><P + Pm) = -1 + ^p(^ H2 + um) ( 64 ) 

Using = UfffPtf,, the latter takes the form 

^e// = -l + ^(^ 2 + ^ 2 -n0)) (65) 
Substituting (45) and (47) into this equation, we obtain 

^'- l+ ^ H2 ~2m + W) + l { 7§ m (66) 

It is now possible to find evolution of the equation of state parameter u e ff(z) with the same 
procedure used in the section 3.1. In fig. lb and fig.2b, the resulting u e ff(z) is plotted for the 
models (40) and (42). Both figures indicate that there are regions in the parameters spaces for 
which crossing the phantom barrier is allowed. 



4 Conclusion 

We have studied phantom behavior for some f(R) gravity models in both Jordan and Einstein 
conformal frames. In Jordan frame, we have used the reconstruction function q(z) fitting to 

^The case of Q — corresponds to ACDM. In this case, the equation (59) states that p m always decreases 
with expansion of the universe. 
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the Gold data set to find evolution of equation of state parameter for some cosmologically 
viable f(R) models. Both models (40) and (42) indicate phantom behavior in some regions of 
the parameters space. 

In Einstein frame, the scalar partner of the metric tensor is separated. Comparing this scalar 
field with the usual notion of a phantom field, we have made an observation that the former 
appears as a quintessence with a minimal coupling to gravity. However, this quintessence field 
interacts with matter sector which allows energy transfer between the two components. This 
interaction plays a key role in the cosmo logical behavior of a particular f(R) model in Einstein 
conformal frame. Then we have used the same reconstruction function for q(z) to study the 
evolution of the effective equation of state parameter. We have found that both models (40) 
and (42) can display crossing the phantom boundary in the same regions of the parameters 
spaces that is used in Jordan conformal frame. 

To compare the behavior of u e ff(z) for each model in the two conformal frames, this function 
is plotted for a common set of the parameters. From fig.l and fig.2, one can see that overall 
behaviors of the equation of state parameters are very similar. In particular, both models 
display crossing the phantom boundary nearly at the same epoch in both conformal frames. 
Despite this similarity, details of the equation of state parameters in the two frames are not 
exactly the same. This result stands against the idea (which is partly presented in [23]) that 
considers both conformal frames with the same physical status. 
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Figures : 



(a) (b) 




Figure 1: The plot of u e ff(z) for the model (40) in a) Jordan and b) Einstein frames. The 
function u e ff(z) is plotted for a common set of the parameters A and n. The figures show that 
crossing the PDL can be realized in both conformal frames. 



(a) (b) 




Figure 2: The plot of u e ff(z) for the model (42) in a) Jordan and b) Einstein frames. The 
function u e ff(z) is plotted for a common set of the parameters 7 and m. The figures show 
that the PDL can be crossed in both conformal frames. 
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